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(Dated: February 14, 2007) 

We present a measurement of the shape of the boson rapidity distribution for pp — * Z/"/* — > 
e + e~ + X events at a center-of-mass energy of 1.96 TeV. The measurement is made for events with 
electron-positron mass 71 < M ee < 111 GeV and uses 0.4 fb _1 of data collected at the Fermilab 
Tevatron collider with the DO detector. This measurement significantly reduces the uncertainties on 
the rapidity distribution in the forward region compared with previous measurements. Predictions 
of NNLO QCD are found to agree well with the data over the full rapidity range. 



PACS numbers: 13.60.Hb, 13.38.Dg, 13.85.Qk 
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I. INTRODUCTION 

Kinematic distributions of Z/-f* bosons produced in 
hadronic collisions provide a wealth of information on 
the fundamental interactions involved. At leading order, 
Z/-f* bosons are produced through the annihilation of a 
quark and an anti-quark, with the partons in the proton 
and anti-proton carrying momentum fractions x\ and X2, 
respectively. In turn, the rapidity of the boson, defined as 
y = 5 In , where E is the energy of the boson and 
Pl is the component of its momentum along the beam 
direction, is directly related to the momentum fractions 

by 

Mz/r ±y 

Here, M z / 7 » is the mass of the boson, and y/s is the cen- 
ter of mass energy. These kinematic distributions can be 
precisely reconstructed when the boson decays leptoni- 
cally since the leptons can be accurately reconstructed, 
and the backgrounds to di-lepton final states are small. 
For low rapidity bosons, the leptons also have small pseu- 
dorapidity, 

j] = - In (tan (6/2)) , 

where 6 is the polar angle and is measured relative to 
the proton beam. High rapidity bosons are more likely 
to have initial states that have maximal \x\ — xi\ for the 
incident partons. 

Although calculations are available at next-to-next- 
to-leading-order in QCD (NNLO) for der/dy for pp — > 
Z/-f* — ► 11 + X [1], few measurements of the differen- 
tial cross section exist [2]. The forward rapidity region 
(\y\ > 1.5) suffers from a smaller cross section and lower 
acceptance than the central rapidity region (|y| < 1.5), 
and has not yet been well tested. The forward region 
probes quarks with low x and high 4-momentum transfer 
squared Q 2 (Q 2 « M§) as well as quarks with very large 
x. Parton distribution functions (PDFs) in this regime 
are mainly determined by jet cross section data, which 
have very different experimental and theoretical system- 
atic uncertainties than the electron measurements pre- 
sented here, and by inclusive lepton scattering data taken 
mostly at much lower Q 2 , which must be evolved to high 
momentum transfer scales using the DGLAP equations 
[3]. 

We measure the normalized differential cross section 



-(-) 



(e x A) avg n; 



obs 



JV total 



where the index i indicates the boson rapidity bin. In the 
first term on the right hand side, e avg is the average effi- 
ciency and A avg is the average acceptance for kinematic 
and geometric cuts. N°^ al is the total number of candi- 
date bosons, and A to ^ al is the total number of background 
events in the sample. In the second term, Aj, N° hs , 



and N^ s are the same as before, but determined in each 
bin i. Aj is the bin width. Dividing by the total number 
of events reduces many of the systematic uncertainties, 
particularly those due to luminosity. 

The DO detector [4] allows efficient detection of elec- 
trons [5] at the large pseudorapidities needed to study 
high rapidity Z/j* bosons. It has a central tracking 
system, consisting of a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT), both located within 
a 2 T superconducting solenoidal magnet, with designs 
optimized for tracking and vertexing [6] at pseudorapidi- 
ties \tid\ < 3 and \r)]j\ < 2.5, respectively. The quantity 
i]D is rj measured from the center of the detector. A 
liquid-argon and uranium calorimeter allows reconstruc- 
tion of electrons, photons, jets, and missing transverse 
energy. The calorimeter is divided into a central sec- 
tion (CC) covering \t]d\ < 1-1 , and two end calorimeters 
(EC) that extend coverage to |?td| « 4.2. Each calorime- 
ter is housed in a separate cryostat [7]. An outer muon 
system, covering < 2, consists of a layer of track- 
ing detectors and scintillation trigger counters in front of 
1.8 T iron toroids, followed by two similar layers after 
the toroids [8]. The luminosity is measured using plas- 
tic scintillator arrays placed in front of the EC cryostats. 
The trigger and data acquisition systems are designed to 
accommodate the high luminosities of Run II. 



II. EVENT SELECTION 

This measurement utilizes a data set of 0.4 fb _1 
collected at the Fermilab Tevatron between 2002 and 
2004. The data are from pp collisions at \/s = 1.96 
TeV. We consider candidate Z/^f* events that decay into 
an electron-positron pair with a reconstructed invariant 
mass 71 < M ee < 111 GcV. The range used is ±20 GeV 
about the mass of the Z boson. 

To optimize the acceptance for electrons at large ry, two 
strategies are used. The first is to require only one of the 
electrons be matched to a reconstructed track. Requir- 
ing a track-matched electron helps to reduce background 
from jets misidentified as electrons, while removing the 
track requirement on the second electron extends the r\ 
coverage beyond that of the tracking system. The sec- 
ond strategy takes advantage of the length of the bunches 
containing the incident protons and antiprotons, which 
has a design length of 37 cm. For our dataset, the z co- 
ordinate [9] of the primary interactions have a Gaussian 
distribution with an rms ranging from 29 ±2 cm to 24 ± 1 
cm. The rms varies with respect to run conditions and 
time. At large values of vertex \z\, some of the decay 
products will travel back through the detector towards 
smaller \z\ values. These particles pass through much 
of the active volume of the tracking system. Typically, 
these events have low background and have the highest 
boson rapidities. 

Events are considered only if a single electron trigger 
fired. The efficiency is (99.0 ± 0.3)% per electron for 
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particles with transverse momentum px > 30 GeV and 
\t]d\ < 2. From these events, di-electron candidates are 
selected by requiring two isolated electromagnetic objects 
which have shower shapes consistent with those of elec- 
trons. A candidate electron is considered to be isolated 
when greater than 85% of its energy is contained within a 
cone of 0.2 in -q — <fi space. Also, they must have deposited 
more than 90% of their energy in the electromagnetic por- 
tion of the calorimeter. One electron must have px > 15 
GeV while the other electron hasp^ > 25 GeV. Electrons 
are defined to be in the CC (EC) region of the calorime- 
ter if they are within \rj D \ < 0.9 (1.5 < \t]d\ < 3.2). In 
the CC region, electrons are not used if they pass near 
EM calorimeter module boundaries. As mentioned pre- 
viously, at least one of the electrons in each di-electron 
pair must be spatially matched to a reconstructed track. 
In addition, all CC-region electrons are required to have 
a track match, so that both legs of a CC-CC event have a 
track match. A total of 19,306 events pass these selection 
criteria. 



III. EFFICIENCIES AND BACKGROUNDS 

Single electron efficiencies are measured from this data 
sample using a "tag-and-probe method." This method 
involves selecting a sample of Z — > e + e~ candidate events 
by applying tight selection criteria to one of the electron 
candidates, the "tag leg," and very loose selection cri- 
teria to the other electron candidate, "the probe leg." 
The reconstructed mass of the tag and probe pair are re- 
quired to be close to that of a Z boson. The tag leg has 
tighter cuts to reduce the amount of background and to 
increase the probability that the event is really a boson 
decay event, and not from jets that are misidentffied as 
electrons. The probe leg has looser cuts and is used to 
test the selection cut in question. While the efficiencies 
are measured with data, Monte Carlo events are used to 
test for biases in the efficiency measurements. For this 
purpose, Z/^j* Monte Carlo events are generated with 
pythia [10] and are processed with a full DO detector 
simulation based on the GEANT software package [11], 
which models the interactions of particles with matter. 
Efficiencies are measured for identification of particles 
like photons and electrons that shower in the electromag- 
netic calorimeter ("EM particles"), shower shape cuts, 
trigger, and track-matching probability. All efficiencies 
are studied as a function of the r\u of the probe electron. 
In addition, some of the efficiencies are measured with 
respect to additional quantities such as px of the probe, 
vertex z position of the event, boson y, or run number. 
Single electron efficiencies are relatively flat in r\x> for the 
CC region, and the values are typically larger than 90%. 
In the EC region, the efficiencies are sensitive to changes 
in the calorimeter geometry, to the finite coverage of the 
tracking system, and to the shape of the distribution of 
event vertices. Due to this sensitivity, effects of variations 
in the width of the vertex z distribution of the course of 



Run II are taken into account with a width-dependent 
efficiency. 

For the acceptance determination, wc use the resbos 
Monte Carlo event generator [12] with CTEQ6.1M input 
PDFs [13, 14]. RESBOS computes the differential cross 
section including NLO QCD corrections and uses resum- 
mation for the low px region. The simulated events then 
are processed with PHOTOS [15] to account for QED fi- 
nal state radiation (FSR). The events then are passed 
through a parameterized detector simulation which has 
been tuned to our data set. To properly apply efficiencies 
in the Monte Carlo, events are weighted based on the rel- 
ative integrated luminosity per run. Figure 1 compares 
data to the Monte Carlo simulation; the simulated signal 
plus background reproduces the data well. The ex A per 
rapidity bin is summarized in Table I. 

The main source of background arises from events with 
jets where one or more of the jets are misidentified as an 
electron. The size of the background is less than 0.8% 
for events where both electrons are detected in the cen- 
tral calorimeter (CC-CC) and less than 6% for the re- 
mainder of the data set. The background for CC-CC 
events is significantly smaller due to the track-match re- 
quirement on both electrons. The jet background is sub- 
tracted by fitting the di-lepton mass distribution with the 
sum of background and signal shapes. The signal shapes 
are taken from the same tuned Monte Carlo as used for 
the acceptance. Two different methods to determine the 
background shape are used. For |y| < 2, the background 
is determined separately for each rapidity bin. The to- 
tal background is measured in the low statistics regions 
of y > 2 and y < — 2 separately. The small numbers of 
events in these areas do not permit the background fits 
to be performed on a bin- by-bin basis. In each of the 
high rapidity regions the background fraction is assumed 
to be constant for subsets of CC-CC, CC-EC, and EC- 
EC events. The background per bin then is determined 
using the number of candidate events per bin collected in 
each subset and the background fraction for that subset. 

Additional background contributions could come from 
events that produce two EM objects in the final state. 
We consider diboson events containing a W plus a W, Z, 
or 7; tt; and Z/j* — > tt events where each r decays to 
an electron. The combined contribution from these addi- 
tional sources is negligible compared to the background 
from jets. The total number of background events per 
bin is presented in Table I. 



IV. SYSTEMATIC UNCERTAINTIES 

A number of contributions to the systematic uncer- 
tainty are considered. These include contributions that 
arise from the determination of the e x A correction and 
from the measurement of the background. 

For estimating the background systematic uncertainty, 
two different background shapes are used in the fits. One 
is obtained from electron-positron events that fail the 
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FIG. 1: Comparisons of data and Monte Carlo plus background are presented for (a) the vertex z distribution and (b) the 
electron-positron invariant mass spectrum. The vertex z plot shows data after all selection cuts. The data in the mass plot 
pass all selection criteria except for the mass cut. Uncertainties shown on the data points are statistical. 



shower shape cuts; the other parameterizes the back- 
ground shape as an exponential curve and incorporates it 
directly into the fit. The exponential fits result in about 
13% more background. The average of the two methods 
is used as the background central value and the differ- 
ence is split and assigned as a systematic uncertainty. 
An additional contribution to the background system- 
atic uncertainty is derived by varying the constraints on 
the signal amplitude used in the background fits and re- 
determining the background. 

The uncertainty on the differential cross section from 
the uncertainties on background ranges from 1.5 — 2.0% 
for \y\ > 0.8. In this region the data mainly come from 
the CC-EC and EC-EC configurations which, because 
only one track-match is required, tends to allow more 
background. For \y\ < 0.8, which is dominated by CC- 
CC data, the uncertainty due to the background is less 
than 1%. 

Several contributions due to the e x A measurement 
are taken into account. These include the uncertainties 
on single electron efficiencies, the electron energy scale 
and energy resolution, the PDFs, and the model of the 
vertex z distribution. 

For single electron efficiencies, several aspects con- 
tribute to the systematic uncertainty. The first two con- 
tributions are derived from data while the third is ob- 
tained using Monte Carlo events. Since efficiencies are 
measured using data, the size of our Z/j* sample in- 
herently has a limited precision. This statistical uncer- 
tainty is included as part of the systematic uncertainty. 
The next component of the efficiencies' systematic un- 
certainty comes from the background subtraction. To 
estimate this contribution, selection cuts are tightened 
on the tag electron to reduce the background at the ex- 



pense of statistical precision. A comparison of the effi- 
ciencies with nominal and tighter cuts is used to estimate 
the systematic uncertainty from the background subtrac- 
tion technique. Lastly, the tag-and-probe method used 
in the efficiency measurement may produce a biased re- 
sult if the efficiency for the probe electron passing the 
selection criteria is correlated with that of the tag elec- 
tron. We estimate the size of this bias with the PYTHIA 
Monte Carlo sample mentioned above, which includes a 
full detector simulation. Efficiencies measured using gen- 
erator level information about the true particle identities 
are compared to the same efficiencies measured via tag- 
and-probe. The difference is used as a contribution to 
the systematic uncertainty. 

For the parameterized detector simulation, the energy 
response and resolutions are tuned using the width and 
peak position of the electron-positron mass distribution 
from the data sample. Kinematic variables that are cor- 
related with the boson rapidity are not used to tune the 
detector simulation parameters. Changing the tuning 
method leads to slight variations in the energy scale and 
resolution parameters. From these variations, we esti- 
mate the contribution due to uncertainties on the elec- 
tron energy scale. The boson rapidity measurement is 
not sensitive to the energy resolution. 

CTEQ6.1M PDFs are defined by twenty orthogonal 
parameters. Each parameter has an uncertainty which 
is shifted separately in the positive and negative direc- 
tion to provide a set of forty PDFs for determination of 
the uncertainty. The acceptance is reevaluated with each 
PDF. Following the prescription presented in Ref. [14], 
we compare each acceptance to that obtained with the 
nominal PDF set. The differences are combined into a 
PDF uncertainty, with a distinction made for sets that 
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FIG. 2: Relative uncertainties plotted against the boson ra- 
pidity. 



FIG. 3: D0 Run II measurement of ^da/dy vs y. The inner 
(outer) error bars show the statistical (total) uncertainty. 



increase or decrease the acceptance. 

As mentioned above, the shape of the vertex z distri- 
bution varies with time. The width of the distribution 
can depend on a number of factors. The beam tuning has 
changed this width over the course of Run II. Also, the 
time elapsed since beam injection can affect the width. 
Since the probability of an electron to have a track- match 
depends in part on the z position of the primary vertex, 
knowledge of the vertex distribution can directly effect 
the acceptance correction. Samples of vertex z distri- 
butions extracted from the data in blocks corresponding 
to different instantaneous luminosities are used to model 
the vertex distribution in the Monte Carlo simulation. 
Selection criteria that produce the widest and narrow- 
est vertex distribution widths are used to estimate the 
systematic uncertainty. 

The main contributions to the total systematic uncer- 
tainty depend on the boson rapidity. At small values of 
|y|, the main sources are the single electron efficiencies 
(« 2%) and the vertex z distribution (< 1%). For mid- 
range \y\, the largest contributions are due to the electron 
efficiencies and the background subtraction. The size of 
each is roughly 2%. In the high rapidity region, \y\ > 2, 
the main sources are from the electron efficiencies, the 
background, and the PDFs. The combined uncertainty in 
this region ranges from 3% to 10% and increases with \y\. 
The relative total systematic uncertainty along with the 
contributions to the uncertainties from the background 
and the e x A are presented in Fig. 2. Contributions to 
the ex A uncertainty also are presented for each rapidity 
bin in Table II. 

To cross check our result, we split the data into in- 
dependent sets based on criteria that should not affect 
the result. These include dividing the data based upon 



(a) time period for data collection, (b) different ranges of 
instantaneous luminosity, and (c) the calorimeter region 
in which the electrons are detected. Cross sections from 
independent subsets are compared to look for inconsis- 
tencies. Subsets in (a) are sensitive to hardware changes 
over the course of the data set and/or changes to the trig- 
ger menu used in collecting the data. Subsets in (b) have 
different vertex z distributions and will not agree if the 
vertex distribution is modelled poorly. Subsets in group 
(c) compare data from three separate calorimeters. All 
of the cross checks give results that are consistent within 
uncertainties. 



V. RESULTS 

A plot of ^der/dy is given in Fig. 3 for Z/j* events 
within a mass range of 71 < M ee < 111 GeV. The inner 
(outer) error bars show the statistical (total) uncertainty. 
In Fig. 4 the result is shown vs \y\. For bin centering we 
follow the prescription given in Ref. [16]. The center of 
the bin is located at the average value of the expected 
distribution over the bin. For this purpose we use the 
NNLO calculation generated from code made available 
from Ref. [1]. This is a small effect and for the two dec- 
imal places of precision used here, the procedure gives 
the same result as using the bin center. Due to the finite 
resolution of the DO detector, some fraction of the events 
in a given rapidity bin originates from a neighboring bin. 
For this analysis, about 5% of the events migrate to each 
adjacent bin. Even though the effect is small, a migra- 
tion correction is included in the e x A determination. 
The uncertainties in -da/dy are dominantly statistical 
for all measured rapidity bins. 



DO, 0.4 fb 1 




Boson Rapidity, |y| 

FIG. 4: D0 Run II measurement of ^da/dy vs \y\. The inner 
(outer) error bars show the statistical (total) uncertainty. The 
curve is a NNLO calculation from [1] using MRST 2004 PDFs. 

The values for the fraction of the cross section in each 
rapidity bin also are listed in Table I. No information on 
the bin-to-bin correlations is included in the table. Since 
the systematic uncertainty is small compared to the sta- 
tistical uncertainty a correlation matrix is not included. 
The curve in Fig. 4 is a NNLO calculation from Ref. [1] 
generated with MRST 2004 NNLO PDFs [17]. The cal- 
culation agrees well with our data, with a x 2 /d.o.f. of 
20.0/27. Our result improves upon previous measure- 
ments over the full range in y, especially in the forward 
region. Figure 5 shows the relative uncertainties from the 
most recent boson rapidity measurement [2] and from this 
analysis plotted vs \y\. For comparison purposes, Fig. 5 
also includes the relative uncertainty due to PDFs for a 
NLO calculation of the differential cross section. This 
curve uses CTEQ6M uncertainty PDFs and code from 
Ref. [1] . The method is the same as that used for the de- 
termination of the PDF uncertainty on the e x A which 
was discussed earlier. 

In summary, we have presented a measurement of 
■^da/dy for Z/j* measured with electron-positron events 
in the mass range 71 < M ee < 111 GeV. The mea- 
surement is the most precise measurement to date. It 
improves upon previous measurements and gives a sig- 
nificantly more precise measurement of the boson rapid- 
ity distribution in the high rapidity region which probes 



the small x, high Q 2 portion of the parton distribution 
functions which is least constrained by other data. The 
fractional uncertainty in the highest rapidity bin is re- 
duced by 30%. We find the result to be consistent with 
a recent NNLO calculation. The current measurement 
is performed with (10-20)% of the expected Run II in- 
tegrated luminosity. An improved result is foreseen with 
DO, 0.4 fb 1 

£• F » I 

c - Relative Uncertainty 

"5 
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FIG. 5: Relative uncertainties from this measurement and 
from the CDF Run I result. Also shown is the PDF uncer- 
tainty on the differential cross section using the CTEQ6M 
uncertainty PDF sets. The values for the CTEQ6M curve are 
generated with code from Ref. [1]. 



the inclusion of additional data, which will reduce the 
current still-dominant statistical uncertainty. 
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TABLE I: Summary of the measurement of per rapidity bin for Z/y* — ► e + e events with mass 71 < M ee < 111 GeV. 
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TABLE II: Contributions to the 
described in the text. 
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